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An Homogeneous Bright QSO Survey (HBQS) has been carried out in the framework 
of an ESO Key programme. 327 QSOs (with Mb < —23 and 0.3 < z < 2.2) have been 



^ \ selected over an area of 555 deg in the magnitude range 15 < i? < 18.75. For magnitudes 

^ ' brighter than B = 16.4 the QSO surface density turns out to be a factor 2.2 higher than 

what measured by the PG survey, corresponding to a surface density of 0.013t;oo5 deg~'^. 
^ . If the data from the Edinburgh QSO Survey are included, an overdensity of a factor 2.7 

^ I is observed, corresponding to a surface density of 0.016 ±0.005 deg~^. In order to derive 

the QSO optical luminosity function we used Monte Carlo simulations that take into 
account of the selection criteria, photometric errors and QSO spectral slope distribution. 
We have combined our data with the Edinburgh QSO Survey, the Large Bright QSO 
Survey, Durham/AAT survey, ESO/AAT faint survey and the (ZM)^B survey. 

The luminosity function can be represented with a Pure Luminosity Evolution, 
L{z) oc (1 + z)'', of a two power-law both for go = 0.5 and go = 0.1. For go = 0.5 
the k evolution parameter is k = 3.26 ± 0.07, slower than the previous Boyle's (1992) 
estimate k = 3.45. A flatter slope /? = —3.72 ± 0.13 of the bright part of the luminosity 
function is also required. If a spread in the QSO spectral slope of cr-y = 0.3 and 0.5 is 
taken into account, the k parameter drops to 3.2 and 3.0 respectively. The observed over- 
density of bright QSOs is concentrated at redshifts lower than 0.6. It results that in the 
range 0.3 < z < 0.6 the luminosity function is flatter than observed at higher redshifts. 



^ Based on data of the ESO Key-Programme "A Homogeneous Bright QSO Survey" 
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In this rcdshift range, for Mb < —25, 32 QSOs are observed instead of 19 expected 
from our best-fit PLE model. This feature requires a luminosity dependent luminosity 
evolution in order to adequately represent the data in the whole 0.3 < z < 2.2 redshift 
interval. The observed overdensity of low redshift bright QSOs could be originated by 
the contribution of a slower evolving population of radio-loud QSOs. 

Subject headings: 159 - Quasi-stellar Objects 
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1. INTRODUCTION 

The study of the evolution of the QSO Luminosity 
Function (LF) is of primary importance in order to 
understand the physic driving the central "engine" 
of AGNs (e.g. Efstathiou & Rees 1988, Haehnelt 
& Rees 1993). Studies of the large scale structure 
of the Universe (e.g. Andreani & Cristiani 1993), 
and of the extragalactic X-ray and UV background 
(e.g. Madau, Ghisellini & Fabian 1994; Comastri et 
al. 1995; Haardt & Madau 1996) are based on the 
knowledge of the optical QSO LF and its relationship 
with other wavebands. The LF of the optically se- 
lected QSOs is also the starting point to disentangle 
the evolutionary properties of the two populations of 
the radio-loud and radio-quiet QSOs (La Franca et al. 
1994). 

The general behaviour of the LF is fairly well es- 
tablished in the redshift interval 0.3 < z < 2.2 
for which color techniques provide reliable selection 
methods (see Hartwick and Schade (1990), Boyle 
(1992) and Hewett and Foltz (1994) for a review of 
the subject). The prevailing model was proposed by 
Boyle, Shanks & Peterson (1988) and was based on 
a sample of 420 faint [h < 20.9) AGNs, combined 
with the Palomar Green (PG) Bright QSO Sample 
(BQS, Schmidt & Green 1983) consisting of 92 QSOs 
with B < 16.4. The model consists of a double- 
power-law shape invariant with redshift, with a steep 
bright end dN/dL oc L~^-^, and a flatter faint end 
dN/dL oc L~^-^. The evolution is simply represented 
by the so-called Pure Luminosity Evolution (PLE) in 
which the characteristic luminosity {L*) of the break 
between the two power-laws evolves as a function of 
redshift according to L*{z) oc (1 -|- z)^-'^. 

In the last years a number of papers have been 
published suggesting that the standard two power-law 
model in the framework of a pure luminosity evolution 
may not be an adequate representation of the data 
(Goldschmidt et al. 1992; Miller et al. 1993; Hewett, 
Foltz & Chaffee 1993; Hawkins & Veron 1995). The 
most controversial points are the amount of flattening 
of the faint part of the LF, first discovered by Boyle 
et al. (1988), and the shape of the bright part which, 
in the last decade, has been based almost entirely on 
the PG bright QSO survey data. 

In this paper we present a new determination of 
the QSO LF and its evolution. We devote special care 
to the shape of the bright part of the LF, for which 
we have derived new data from a bright QSO sam- 



ple: the Homogeneous Bright QSO Survey (HBQS, 
partly published in Cristiani et al. 1995). In the 
computation we take also into account the most rele- 
vant already existing QSO samples, in particular two 
new bright QSO surveys: the Edinburgh QSO Sur- 
vey (EQS, Goldschmidt et al. 1992), and the Large 
Bright QSO Survey (LBQS; Hewett, Foltz & Chaffee 
1995 and references therein). 

2. THE SAMPLES 

As discussed in section 4, we have used Monte 
Carlo simulations in order to derive the best fit QSO 
LF. For this reason our analysis is mainly based on 
a set of QSO samples for which the selection criteria 
and the photometric properties are well defined or at 
least can be acceptably reproduced. In addition to 
the knowledge of the K-corrcction and QSO colors, 
for a rigorous determination of the LF we need the 
definition for each sample of: 

• the flux limits and selection criteria; 

• the distribution of the photometric errors as a 
function of apparent magnitudes; 

• the galactic absorption; 

• the color equations of the photometric system. 

The QSO LF will be computed in the Johnson B 
absolute magnitudes. In the following color equations 
the "natural" photometric systems are marked with 
primed letters, while the Johnson/Cousins system is 
represented with unprimed letters. 

2.1. The HBQS 

The Homogeneous Bright QSO Survey (HBQS) 
was started in 1989 in the framework of an ESQ Key- 
programme. The survey covers a total of 555 deg^ 
subdivided in 22 ESO/Schmidt or UKST fields at 
high galactic latitude around the south galactic pole 
{b < -60°). 

Two Schmidt plates for each bandpass U (IIa-0 
+ UGl or Illa-J + UGl), B' (IIa-0 + GG385) or 
Bj (Ilia- J + GG395), V (Ila-D + GG495), R (Illa- 
F + RG630) or OR (IIIa-F + OG590), and / (IV- 
N -I- RG715) have been obtained at the UKST or 
ESO/Schmidt telescope, usually within a few months 
interval, in order to minimize the effects of variability. 
The plate material has been scanned on the COSMOS 
microdensitometer. The resulting tables, containing 
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the instrumental magnitndcs and other useful param- 
eters for the objects detected in each plate, have been 
merged together in one table per field. Only ob- 
jects with at least 4 detections (in 10 plates) have 
been included in this final table. According to the 
ESO/UKSTU numeration, the 13 ESO fields which 
have been selected in the B' photometric system are: 
290, 349, 355, 406, 407, 408, 410, 468, 469, 470, 474, 
479, 534. The 9 UKST fields which have been selected 
in the Bj photometric system are: 287, 295, 296, 297, 
351, 411, 413, SA94, SGP. 

We have selected as candidates all the UVx "not- 
extremely-extended" objects satisfying a type of mod- 
ified Braccesi less-restricted two-color criterion (La 
Franca, Cristiani & Barbieri 1992). The principal 
color for the selection has been choscin to be the U—B' 
or the U — Bj according to the type of the available 
plate material. The secondary color has been chosen 
to be Bj — i? or Bj — V' or B' — V' , always preferring 
the combination giving the most reliable photometric 
accuracy. 

The adopted color equations (Blair & Gilmore 
1982) are, for the B' (IIa-0 + GG385) system: 

B = B' + 0.11{B-V); 

while for the B'j (Illa-J -\- GG395) system: 

B = B'j + 0.28{B -V) (1) 

The magnitude intervals in which the selection is 
virtually complete vary from field to field between 
15.0 < B < 17.3 and 15.0 < B < 18.8. The 
sample extracted from these data includes 327 QSOs 
(Mb < -23, qo = 0.5, and Hq = 50 Km/s/Mpc). 
Cristiani et al. (1995) describe in detail the pro- 
cedure followed in the data acquisition and reduc- 
tion and report the first results obtained in six fields 
(^ 150 deg^, the ones reaching fainter limiting mag- 
nitudes), providing significant statistics in the range 
(16.75 <B < 18.75). 

2.2. The Edinburgh QSO Survey (EQS) 

The Edinburgh bright UVx QSO Survey (EQS) has 
been carried out over 13 UKST fields (the area is 333 
deg^ large) with the UBjVRI multicolor system (see 
Goldschmidt et al. 1992). The data reduction was 
very similar to that of the HBQS. The photometric 
rms errors are of 0.09 mag in each bandpass in the 
magnitude interval 15 < i? < 18 were most of the 



QSOs are. The selection criteria is fairly well repre- 
sented by the equations: 

U -Bj < -0.35 and Bj-R>0. 

The intervals of magnitude for which a reliable 
sample can be extracted vary from 15.0 < B < 17.5 to 
15.0 < B < 18.5 in the different fields (Goldschmidt 
1994). For our computation of the LF we have used 
the bright {B < 16.5) subsample of 8 QSOs with 
z > 0.3 published by Goldschmidt et al. (1992). In 
this subsample the AGN with h = 15.79 and z = 0.380 
has been excluded because of a wrong redshift iden- 
tification (Goldschmidt 1994). 

2.3. The LBQS 

The Large Bright QSO Survey covers an effective 
area of 454 deg^ (Hewett, Foltz and Chaffee 1995 and 
references therein) . The photometry has been carried 
out in the Bj bandpass and the errors are on aver- 
age of 0.10 magnitudes, taking also into account zero 
point errors and variations across UKST fields. The 
sample is based on 1055 QSOs and AGNs, with ap- 
parent magnitudes in the range 16.5 < Bj < 18.85 
and redshifts in the range 0.2 < z < 3.4. The se- 
lection has been carried using objective prism plates, 
looking for emission line and/or blue continuum ob- 
jects. Because of the selection of the blue continuum 
objects the sample is complete in the redshift range 
0.2 < 2; < 2.2. It is not possible to directly simulate 
the selection criteria of this survey, however in order 
to check the effects in the LF by the inclusion of such 
a large database, we have temptatively simulated the 
blue continuum selection as an UVx selection: 

U-Bj < -0.30. 

As for the EQS and HBQS samples the adopted color 
equation is (1). In our computation of the LF we have 
restricted the LBQS to a subsample of 220 QSOs in 
the magnitude range 16.5 < B < 18.0. 

2.4. The SA94 Survey 

The SA94 survey is an UVx (UBV) selected sample 
over an area of 10 deg^. The QSOs have been selected 
according to the less restricted Braccesi UVx criterion 
(La Franca aZ. 1992). The completeness magnitude 
range is 15.0 < B < 19.9. The complete sample ex- 
tracted from the data includes 97 QSOs. The B mag- 
nitudes are directly published in the Johnson/Cousins 
system and no color transformation is needed. 
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2.5. The Durham/ AAT Survey 



This sample comes from an UVx stellar object se- 
lection from UK Schmidt U and Bj plates. The Bj 
plates have been calibrated on the B magnitude sys- 
tem. Therefore the magnitudes of this surveys are 
referred as b magnitudes for which the color equation 
is (Boyle et al. 1990): 

5 = 6 + 0.23(5-^-0.9). 

The whole catalogue containing 420 AGNs is sub- 
divided in 34 circular fields having a diameter of 40', 
corresponding to a total area of 11.9 deg^. The com- 
pleteness magnitude ranges vary from field to field, 
the faintest limit magnitude reaching b = 20.9. Dif- 
ferent UVx selection limits have also been adopted in 
the various fields (see Boyle et al. 1990, Table 2). On 
average they correspond tou — b< —0.40. In the red- 
shift interval 0.3 < < 2.2 the sample consists of 353 
QSOs. 

2.6. The ESO/AAT 

This is a faint UJRI survey of 66 AGNs in the 
fields 855, 861, 864 (Boyle, Jones and Shanks 1991), 
where J' is a Ilia- J + GG385 system. The color equa- 
tion is (Kron 1980): 

B = J' + 0.23{B-V). 

The survey covers an effective area of 0.85 deg^ from 
J = 18.0 down to J = 21.75-22.0, and the complete- 
ness redshift range is 0.6 < z < 2.9. The typical rms 
scatter in the calibration is 0.09 mag for the J band- 
pass, 0.15 mag for R passband, and 0.12 mag in the 
U bandpass. The candidates were selected among all 
the stellar objects lying blueward of the stellar locus 
in either U — J oy J — R diagrams. At z < 2.2 the 
subsample consists of 41 QSOs whose selection is well 
represented with 

U -J < -0.30. 

2.7. The (ZM)2b Survey 

This is a faint UVX survey of 54 AGNs from UJF 
plates taken at the ESO 3.6m telescope (Zitelli et al. 
1992). The survey cover an area of 0.69 deg^ down 
to J = 20.85 and 0.35 deg^ down to J = 22.0. In 
the redshift range 0.3 < z < 2.2 the complete sample 
consists of 19 QSOs. The selection criterion is well 
represented by 

U-J< -0.30. 




a 
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Fig. 1.— The QSO number counts for 0.3 < z < 2.2. 
The continuous line represents the best fit LDLE 
model (model C in Table 2), the dashed line rep- 
resents our best fit PLE model (model B in Table 
2), while the dotted line represents the best fit PLE 
model (model A in Table 2) by Boyle et al. (1992). 
The magnitudes have all been corrected for galactic 
extinction according to Burstein and Heiles (1982). 

As in the previous survey the color equation is: 
B = J' + 0.23{B-V). 

3. THE QSO COUNTS 

The B number-magnitude counts for all the sam- 
ples used in the computation of the optical LF are 
shown in Fig. 1. In the following we will re- 
strict our analysis to 0.3 < ^; < 2.2 because this 
is the only interval in which the selection criteria 
allow a reliable and homogeneous processing of the 
data. The magnitudes have all been corrected for 
galactic extinction according to Burstein & Heiles 
(1982). The QSO surface densities from the HBQS 
and EQS are also shown in Table 1. As first in- 
dicated by Goldschmidt et al. (1992), for magni- 
tudes brighter than B = 16.4 the QSO surface den- 
sity turns out to be a factor 2.5 higher than what 
measured by the PG survey (Schmidt & Green 1983). 
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The data from the HBQS collect 7 QSOs at mag- 
nitudes brighter than B = 16.4, corresponding to a 
surface density of 0.013l;[]5]^ deg''^. With the ad- 
dition of the EQS, the surface density at _B < 16.4 
becomes 0.016 ± 0.005 deg~^, corresponding to a to- 
tal of 14 QSOs over an area of 888 deg'^. The PG 
sample lists 44 QSOs at _B < 16.4. In this range the 
surveyed area changes rapidly as a function of the 
limiting magnitude. Combining the number of dis- 
covered QSOs with the areas according to Table 1 of 
Green, Schmidt & Liebert (1986), a QSO surface den- 
sity of 0.006 ± 0.001 deg-'^ is obtained, a factor of 2.7 
lower (significant at 2(t level) than the EQS+HBQS 
estimate. The presence of incompleteness in the PG 
survey has been investigated by Wampler and Ponz 
(1985). It probably originates for QSOs with U — B 
colors near the selection limit because of large rms 
errors. This effect would lead to incompleteness that 
varies with redshift, with the consequence that, in 
particular in the range 0.4 < z < 0.7, it is system- 
atically underpopulated at a given magnitude. This 
hypothesis appears corroborated by the fact that the 
PG survey, in the area covered to a limit B = 16.4, 
has discovered 46 QSOs with z < 0.3, in addition to 
the 44 with z > 0.3 that are relevant for the compar- 
ison with the current sample. For these reasons we 
have excluded the PG sample from our subsequent 
analysis of the evolution of the LF. At fainter mag- 
nitudes, in the interval 16.5 < B < 18.5, the HBQS 
shows a fairly good agreement with the counts of the 
LBQS. 

4. THE LUMINOSITY FUNCTION 
4.1. Deterministic Derivations 

The QSO LF $(M, z) is the measure of the number 
of QSOs per unit of commoving volume with magni- 
tudes in the interval {M, M + dM) as a function of 
redshift z. In the redshift interval 0.3 < z < 2.2 the 
general behaviour of the QSO LF has been usually 
parameterized with a double power-law Pure Lumi- 
nosity Evolution (PLE) model <^{Mb,z) = 

$* 

lQOA[MB-MB{z)]{a+l) J. IQ0A[Mb-Mb{z)KI3+1) 



specified by the redshift dependence of the break mag- 
nitude, 

Mb{z) = M*b- 2.5fclog(l + z), 

corresponding to a power-law evolution in the optical 
luminosity, L* oc (1 -|- z^. <^{Mb-iZ) is expressed in 
units of Mpc~^mag~^ (see Boyle et al. 1988). 

The empirical way of deriving the LF is by col- 
lecting all the available QSO samples, computing 
the transformations to a common photometric band- 
pass, and measuring the QSO space densities in Az, 
AM bins using the 1/Va estimator of the "coherent" 
method of Avni & Bahcall (1980); where Va is the ac- 
cessible volume computed for every object over all the 
samples and not just the sample in which the object 
originated. This method has a number of disadvan- 
tages as it is affected by: 

1. the uncertainties due to the binning of the data. 
As the number of QSOs in the available com- 
plete samples is still small, it is necessary to use 
relatively large bin sizes on the Mb — z plane 
within which the evolutionary effects are not 
negligible. In this way the non-linear dependen- 
cies of the QSO LF on Mb and z cause spurious 
trends on the derived shape of the LF. 

2. the influence of the photometric errors. Cav- 

alierc, Giallongo and Vagnetti (1989) pointed 
out how the uncertainties introduced in the sta- 
tistical samples of QSOs by photometric vari- 
ances and from the fluctuation of the sky vari- 
ations cause an overestimation of the evolu- 
tion timescales (see also Giallongo and Vagnetti 
1992, and Francis 1993). 

3. the spread in the assumed mean spectral index 
of the QSOs used to compute the K-correction. 
If the contribution of the emission lines to the 
broadband flux can be neglected, the formula 
relating the absolute and apparent blue magni- 
tudes is (Schmidt & Green 1983) Mb = 

B-5 \og[A{z)c/Ho] + 2.5(1+7) log(l + z) - 25, 



where A[z)c/Ho is the bolometric luminosity 
where a and /? correspond to the faint-end and bright- distance, and the extrapolation from the ob- 

end slopes of the optical LF, respectively. With this served to the rest-frame B band assumes a 

parameterization the evolution of the LF is uniquely power-law representation of the QSO contin- 

uum oc When quantifying the evolu- 

tion of the population between two redshifts, 
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if the extrapolation is large, it is important to 
know the spread of the distribution of 7. Any 
distribution of the mean spectral index is re- 
flected in an uncertainty in the absolute magni- 
tudes, which increases with redshift and intro- 
duces an artificial luminosity evolution of the 
population. The dispersion in the power-law in- 
dex 7 in the rest-frame ultraviolet is substantial 
(Sargent, Steidel and Boksenberg 1989; Francis 
et al. 1991; Schneider et al. 1991), with esti- 
mates of cr-y = 0.3 — 0.6. For a sample of QSOs 
measured at the same observed wavelength the 
differential effect between two redshifts zi, zi is 

AMb = 2.5A7log[(l + ^2)/(l + ^1)]. 

In the redshift range of our analysis 0.3 < ^; < 
2.2, changing from z = 0.3 to z = 2.2 results in a 
uncertainty (Tmb —o.^.K value A7 = 0.5 would 
correspond to an error in the space density of a 
factor 3 in the bright part of the LF assuming 
a power-law index /3 = —3.5. In the case of 
two QSOs of the same redshift and apparent 
magnitude B but different spectral index, the 
difference in absolute magnitude is 

AMs = 2.5A7log(l -h z). 

Changing from 7 = to 7 = 1.5 would produce 
at z = 0.4 a difference in absolute magnitude 
of AMb = 0.4 mag. At 2; = 2.2 the difference 
becomes 1.9 mag corresponding to an error in 
the space density of a factor 100 in the bright 
part of the LF (assuming a power law index 
/3 = —3.5). The final result of this effect is a 
blurring of the LF. 

The "standard" maximum likelihood techniques 
(Marshall et al. 1983, Boyle, Shanks and Peterson 
1988) do not need binning of the data and overcome 
the disadvantage described in item (1), but do not 
solve problems (2) and (3). 

In order to take into account also these latter two 
effects, it has been proposed to compute a detection 
probability as a function of apparent magnitude, red- 
shift and Spectral Energy Distribution (SED) (see 
Warren, Hewett & Osmer 1994). The detection prob- 
ability p(m, z, SED) is derived computing the frac- 
tion of QSOs that would have been selected by each 
survey according to the selection criteria, photometric 
errors, and the spectral index. Thus, for each redshift 



and magnitude interval, the space density is evaluated 
as SI/ J p{z)dVa. Eventually, in order to determine 
the range of parametric luminosity-function models 
that arc consistent to the data, the observed and ex- 
pected number of QSOs are compared through max- 
imum likelihood or least techniques. 

4.2. The Monte Carlo Simulations 

The above described procedures do not lead to a 
fully satisfactory and rigorous determination of the 
LF. The photometric errors and the spread in the 
QSO SED properties make the nature of the prob- 
lem intrinsically stochastic. Even taking into account 
the influence of different SEDs in the computation 
of the absolute luminosity, the problem remains ba- 
sically unsolved. For example, because of the effects 
of the photometric errors, two QSOs observed at the 
same apparent magnitude and with the same SED 
and redshift may be produced by objects of different 
absolute magnitudes. This means that it is impossi- 
ble to establish a bi-unique (one-to-one) relationship 
between the observed apparent magnitudes and the 
absolute magnitudes as a function of z and SED. 

For these reasons we have followed a different ap- 
proach, simulating, with Monte Carlo techniques, 
what happens when a selection of a sample is made. 
For each sample the following procedure has been ap- 
plied: 

1 . Assuming a theoretical LF distribution $ (Mb ,z), 
a number (proportional to the surface area of 
the survey) of QSOs (i.e. pairs Mb,z) have 
been randomly extracted; 

2. The apparent magnitudes of the QSOs have 
been computed following the K-correction ac- 
cording to Cristiani & Vio (1990), and deriving 
the average galactic absorption Ab from the HI 
maps of Burstein and Heiles (1982). The result- 
ing apparent B magnitude is: 

B = MB+5loglA{z)c/Ho]+25+K{z)+AB+e{B) 

where e{B) simulates a noise with Gaussian dis- 
tribution with variance 

O"! = <^lhot{B) + cr|^^. 

(7^ takes into account both the photometric er- 
rors of the survey (Jphot{B) and the apparent 



7 



magnitude dispersion due to the QSO spectral 
slope dispersion a-y 

crs,^, = 2.5cr^log(l + z). 

3. The apparent colors have been generated follow- 
ing the average QSO colors-redshift dependence 
-^1,2(2) (Cristiani & Vio 1990). Between each 
couple of bandpasses centered at Ai and A2 the 
resulting apparent color is: 

Ci,2 = Fi,2iz) + ep,,^ 

where 6^1,3 simulates a noise with Gaussian dis- 
tribution with variance 

2 takes into account the photometric errors 
of the two bandpass tri and (T2 respectively, and 
the QSO color dispersion due to the spread 
in the slope of the spectra 

ap-y = 2.5cr-Y log — ^ . 

Ai 

Note that the error on the color crF,7 is corre- 
lated with the error on the apparent magnitude 
(7B.-I- For each simulated QSO the same value 
of 7 has been used to compute the apparent 
magnitude and colors. 

4. The flux limits and selection criteria of the sur- 
vey have been applied to the resulting mag- 
nitudes and colors in order to select the "ob- 
served" QSOs. 

In addition to statistical errors, the presence of sys- 
tematic errors in the QSO photometry could induce 
biases in the Monte Carlo simulations. In particu- 
lar the applicability of the Blair and Gilmore (1982) 
transformations, which have been calibrated on a 
range of colors (—0.4 < U — B < 2.0), not entirely rel- 
evant for low-redshift QSOs, is not straightforward. 
It is reassuring that the object for which the effect 
might be significant for the selection/non-selection, 
are those closer to the stellar locus, where the color 
transformations are still reliable. 

The observed QSO magnitudes and colors are also 
affected by variability that in principle gives origin to 
at least two effects: 



1. The shape of the LF is modified (blurred) by 
variability. We will derive the QSO LF as it is 
observed at a given instant, so we will neglect 
the effect from this point of view. How such a 
"snapshot" LF is translated into an "average" 
LF, describing the distribution of the average 
luminosities of QSOs on long timescales, is mat- 
ter for future variability studies. 

2. Variability is going to affect the observed colors 
of QSOs, acting as an additional photometric er- 
ror, if the photometric material in the different 
colours is obtained at different times. The criti- 
cal timescale above which the effects of variabil- 
ity become non-negligible can be set at about 
one year in the observer's rest frame (Cristiani 
et al. 1996). However, the effects of variability 
can be modeled with a spectral slope variable 
in time. In particular this can reproduce the 
observed dependence of the variability on the 
wavelength (Di Clemente et al. 1996, Cristiani 
et al. 1996). In this a variable spectral slope 
could manifest itself in an ensemble of QSOs as 
a spread in the observed spectral index, that is 
already included in our computations. There- 
fore we have not introduced further effects of 
the variability on the spread of the QSO colors. 

4.3. Results 

As discussed in section 2, in order to compute the 
best fit LF we have added the QSOs of the HBQS and 
EQS to all the optically selected samples for which it 
was possible a reliable assessment of the selection cri- 
teria and the photometric quality. In this way, a to- 
tal sample of 1022 QSOs has been obtained. The fits 
have been carried out by minimizing the statistics 
derived from the comparison of the obsc;rved {B, z) 
distribution with 2000 simulations of each theoreti- 
cal LF model. Table 2 summarizes the results for a 
number of different models. 

The errors quoted for the parameters are 68 per 
cent (1 cr) confidence intervals. They correspond to 
variations of Ay^ = 1.0, obtained perturbating each 
parameter in turn with respect to its best-fit value, 
and looking for a minimization with the remaining 
parameters free to float. 

The significance of the fitting has also been tested 
for goodness-of-fit in the lowest redshift interval 0.3 < 
z < 0.6, in which the fitting probabilities have been 
computed using the two-dimensional Kolmogorov- 
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Smirnov (2D KS) test by Fasano and Franccschini 
(1988). See Press et al. (1992) for a complete de- 
scription of the algorithm. 

In the following we first focus our discussion on a 
Qo = 0.5 universe, we will later see how our results 
change for Qo =0.1. 

In Fig. 2 our best fit PLE model (B) is shown. The 
data points, representing the estimates of the QSO 
space densities, are compared to the fitted LF. The 
observed QSO space densities have not been directly 
computed through the l/Va technique, that, due to 
the binning of the data, can give origin to spurious 
trends on the shape of the LF (see sub-section 4.1). 
For this reason we have preferred to exploit our Monte 
Carlo simulation code in order to visualize the LF. 
The value of each observational point in Fig. 2 has 
been computed by comparing, in each absolute mag- 
nitude/redshift bin, the number of observed QSOs, 
No, with the number of QSOs simulated with the fit- 
ting LF, iV^. Let us call <l>o and the corresponding 
luminosity functions and \1/(M, z) the selection func- 
tion. In each bin we have: 

No _ J^{M,z)^o{M,z)dMdz 
~ J^{M,z)^s{M,z)dMdz' 

where the integrals are computed over the correspond- 
ing bin. In our case, in each bin, we can assume that 
the shape of the observed and fitting LFs are very 
similar, apart from a constant factor k, which varies 
from bin to bin: k^o- Using this assumption 

we can write: 

N^ ^ ^o{M,z) 

for all values M and z in the bin. Consequently, for 
proper visualization of the data, the observed QSO 
densities have been computed by multiplying the cor- 
responding value of the simulating LF, at the cen- 
tral magnitude and redshift of the bin, by the ratio 
between the number of observed and simulated QSOs 
in the corresponding bin: 

TV 

^o{M,z)c^^^s{M,z). 

The 1 a upper and lower limits have been estimated 
with the Poissonian statistics uncertainties on No- 

As described above, the LF fitting technique is com- 
pletely independent from these assumptions, being 
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Fig. 2. — The QSO luminosity function. The points 
connected with a continuous line represent the obser- 
vations (see text for details about their derivation). 
The dashed lines correspond to the best-fitting Pure 
Luminosity Evolution (model B in Table 2). Error 
bars are based on Poisson statistic and correspond to 
68 per cent (Itr) confidence intervals. 

based on a minimization of the difference between 

No and Ns in bins of B and z. 

Our PLE model B gives a satisfactory fit of the 
{B,z) distribution with a global probability of 
0.21. The PLE model A by Boyle (1992), compared 
with the present data, provides a lower but still ac- 
ceptable probability of 0.12, however with one addi- 
tional parameter, Zcul = 1-9, representing the redshift 
at which the luminosity evolution "switches off" . At 
redshift greater than Zcut no further luminosity evo- 
lution takes place. Our B model has a flatter bright 
slope /? (—3.7 compared to —3.9) than the A model. 
The evolution parameter k is also smaller (3.26 com- 
pared to 3.45). 

This difference in the parameter /3 is originated 
by a larger density of luminous QSOs (expecially at 
z < 0.6) in comparison with the previous data derived 
from the PG sample (compare Fig. 2 with Fig. 3b 
of Boyle 1992). The higher evolutionary rate k deter- 
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mined by Boyle (1992) is a result of the introduction 
of the redshift cut off at Zcut = 1-9, which with our 
data results unnecessary. 

Both A and B models provide an inadequate sim- 
ulation of the distributions of the observed data for 
the low redshift domain, 0.3 < 2; < 0.6. For this 
subsample the 2D KS statistics test rejects model 
A at the 0.01 level and model B at 0.02 level. At 
magnitudes brighter than Mb = —25, in the interval 
0.3 < z < 0.6, 32 QSOs are observed, while 16 and 19 
QSOs are predicted by model A (a 4 ct discrepancy) 
and B (a 3 (7 discrepancy) respectively. But as the low 
redshift {z < 0.6) subsample contain only 5 per cent 
of the complete data set, in a global comparison of the 
whole data sample in the interval 0.3 < z < 2.2, the 
models follow the evolution of the larger fraction of 
QSOs at higher redshift, allowing the probability 
of the global fit to reach satisfactory levels. 

No significant difference is obtained by fitting the 
data excluding each of the three bright samples (HBQS, 
EQS and LBQS) in turn (models D, E and F). 

It is interesting to see how the inclusion of a spread 
in the theoretical average QSO spectral slope modifies 
the best fit luminosity function. As expected (Gial- 
longo and Vagnetti 1992), larger values of the spread 
of the slope correspond to slower luminosity evolution 
parameters and steeper luminosity functions (model 
G with = 0.3, and model H with u-y = 0.5). 

In order to fully represent the data in the redshift 
range 0.3 < z < 2.2 we tried several functional mod- 
ifications of the PLE. As first step we tried, without 
significant results, to add one or at maximum two pa- 
rameters in order to take into account of a dependence 
of the bright slope /? with redshift. 

The best description of the observed data has been 

obtained by decreasing at low redshift the luminosity 
evolution of the bright QSOs. This has been obtained 
by including a dependence on the redshift and abso- 
lute magnitude of the evolution parameter k such as: 

for Mb<M* : k = ki + ^2(^3 - M*)e-^/-(§) 
for Mb> M* : k = ki 

where M* is the magnitude of the break in the two 
power-law shape of the LF. This luminosity depen- 
dent luminosity evolution (LDLE) model (model C in 
Table 1) has resulted in a better fit of the data (see 
Fig. 3) giving a test probability of 0.49 in the 
whole B,z plane, and an acceptable 2D KS proba- 
bility of 0.09 in the redshift interval 0.3 < ^ < 0.6. 
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Fig. 3. — The QSO luminosity function. The points 
connected with a continuous line represent the obser- 
vations (see text for details about their derivation). 
The dashed lines correspond to the best-fitting Lumi- 
nosity Dependent Luminosity Evolution (model C in 
Table 2). Error bars as in Fig. 2. 



With this model, in this redshift interval and for mag- 
nitudes brighter than Mb = —25, 29 QSOs are ex- 
pected in comparison with the 32 observed. 

In this way, as shown in Fig. 1, the LDLE repro- 
duce the higher counts of bright QSOs at B < 16.4, 
discovered by the HBQS and EQS surveys, much 
better than the PLE models: a surface density of 
0.017 deg~^ is predicted with respect to 0.016 deg~^ 
observed. The A model by Boyle (1992) predicts a 
surface density of 0.008 deg~'^. 

It is important to notice that both the PLE and 
the LDLE models do not reproduce a 2tj significant 
observed underdensity of faint QSOs (—23 > Mb > 
—25) in the redshift range 0.3 < 2: < 0.6. In this 
bin, with the PLE B model, 125 QSOs are observed 
and 161 expected. It is possible to find less signifi- 
cant indications of this feature in the data analyzed 
by Boyle (1992). We examined the possibility that 
this observed underdensity of the faint part of the lu- 
minosity function at low redshift is due to an underes- 
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timated incompleteness of the UVx optical selection 
method at low redshift. For this reason we have car- 
ried out simulations assuming a redder U — B QSO 
color of 0.2 magnitudes in the interval 0.3 < z < 0.8. 
The effects were still insufficient to account for the 
underdensity at low redshift (125 QSOs observed and 
154 expected). As our best fit probabilities are sat- 
isfactory, we preferred not to model the low redshift 
faint feature. It is worth to notice that the determina- 
tion of the local (z = 0) AGN (Seyfert 1 plus QSOs) 
LF from the Hamburgh Bright QSO Survey (Wisotzki 
et al. 1996a) shows similar trends, with evidences 
of higher densities at bright magnitudes and lower 
densities at faint magnitudes in comparison with the 
de-evolved PLE LF by Boyle (1992) (Wisotzki et al. 
1996b; Kohler 1996). 

Adopting a Qo = 0.1 Universe, as shown in Table 
2, the fitting of the data with the PLE model (1) by 
Boyle (1992) is rejected at 0.004 confidence level. Our 
PLE model (L) obtains an acceptable global represen- 
tation of the data (x^ probability 0.36), slightly better 
than the PLE fit in a go = 0.5 Universe. However, in 
the 0.3 < z < 0.6 redshift range the PLE still un- 
derestimates the expected number of QSOs, although 
the discrepancy is less serious than in the Qo — 0.5 
case: 44 QSOs with Mb < —25 are observed versus 
33 predicted in our L model (24 in the I model). The 
introduction of a LDLE parameterization according 
to Eq. |2| provides satisfactory representation of the 
observations in the low-z domain also for Qo = 0.1. 

5. DISCUSSION 

What is the origin of the flattening of the lumi- 
nosity function at low-redshift and is it confirmed by 
observations in other wavebands? 

In the last years evidence has been produced that 
optically selected QSOs have a bimodal distribution 
of radio flux densities, and that the fraction of ra- 
dio loud QSOs decreases with increasing redshift and 
decreasing absolute optical luminosity (Miller et al. 
1990; Visnovsky et al. 1992; Schneider et al. 1992; 
Padovani 1993; La Franca et al. 1994). 

La Franca et al. (1994) have determined that the 
radio-loud fraction is substantially higher for QSOs 
with z < 1, a result essentially due to the higher 
overall radio-loud fraction in the predominantly low 
redshift PG sample. The radio-loud fraction has also 
been found to be larger at bright absolute magni- 
tudes, regardless of whether or not the PG sample 



is included in the analysis. La Franca et al. (1994) 
have been able to reproduce the data with a radio- 
loud optical luminosity function (OLF) characterized 
by a lower space density, a similar shape and a smaller 
evolutionary parameter k (2.7 < k < 3.1 ) with re- 
spect to the parent total (radio-quiet plus radio-loud) 
QSO population. According to this picture the flat- 
tening of the bright part of the OLF corresponds to 
an increase of the radio-loud OLF contribution. 

However, the decrement of the radio-loud fraction 
with redshift relies essentially on the PG sample and 
has not been confirmed by a study of 256 QSOs in 
the LBQS (Hooper et al. 1995). Nevertheless there is 
a very small overlap in the magnitude redshift plane 
between the LBQS and the PG sample: the LBQS 
QSOs are too faint and too high-redshift to be com- 
pared with the PGs. Our observed overdensity of low 
redshift (0.3 < z < 0.6) QSOs involves the QSOs with 
B < 16.5 and Mb < -25 but in the LBQS there are 
only two such QSOs with B < 16.5, to be compared 
with the 45 QSOs of the PG sample. In the EQS plus 
HBQS combined sample there are 16 QSO brighter 
than 16.5. It becomes crucial to obtain radio data of 
this QSO subsample. 

Other information on the evolution of the LF and 
on the dependence of the radio loud fraction on mag- 
nitude and redshift can be derived from QSO samples 
selected in the soft X-rays. 

A study of the EMSS AGN sample by Maccacaro 
et al. (1991) and Delia Ceca et al. (1992) has shown 
that a PLE model with a power-law dependence on 
the redshift, Lx{z) = Lx{0){l + z)'^ , provides the sim- 
plest description of data. Their best-fitting evolution 
rate {kx = 2.56) is almost one unit lower than the 
optical best fit. Boyle et al. (1993) used a new sam- 
ple of faint AGNs selected by ROSAT, combined with 
the EMSS sample, to derive a value of A; = 2.8. Using 
the same data Franceschini et al. (1994) showed that 
the evolution rate of X-ray selected AGNs could have 
been underestimated and actually be comparable to 
that of the QSOs in the optical, with a linear scaling 
of the emission in the two bands. Boyle et al. (1994) 
using an enlarged sample of 107 faint AGNs selected 
by ROSAT, combined with EMSS sample, derived an 
evolution L^iz) = Lx{0){l + z)^-^^ at 2 < 1.6, quite 
similar to that of optically selected QSOs. The lat- 
ter two analyses favour a picture in which the same 
QSO population is observed, both in the soft-X and 
in the optical. If this is true, at low z, bright X- 
ray selected samples should show a flattening of the 
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bright part of the OLF and an increase of the ra- 
dio loud fraction. Unfortunately the brightest X- 
ray selected QSO sample, the EMSS, is not bright 
(large) enough to probe this domain of the OLF. In 
the interval 0.1765 < z < 0.4286 the brightest signif- 
icant bin is at LogLx = 45 ergs~^ (0.3 — 3.5 KeV) 
(Maccacaro et al. 1991), which corresponds to about 
Mb = —25 (go = 0.1, see Fig. 5 in Franceschini 
et al. 1994, and La Franca et al. 1995), in the 
interval 0.4286 < z < 08182 the EMSS reaches 
LogLx = 45.5 ergs~^ (0.3 — 3.5 KeV), correspond- 
ing to about Mb = — 26. It is interesting to notice 
that Jones et al. (1996), at bright magnitudes in the 
interval 0.0 < 2; < 0.4, show a flatter X-ray LF with 
respect to their PLE model (see their Figure 5). 

Delia Ceca et al. (1994) have analyzed the radio 
properties of the EMSS: they find that in a subsam- 
ple of optically luminous QSOs, 30% of the objects 
are radio-loud, while in the complementary faint sub- 
sample the fraction of radio-loud drops to 3%. No 
evidences are found that the evolutionary proper- 
ties of the radio-loud objects differ significantly from 
the radio-quiet ones. However, the shape of the de- 
evolved XLF of the two classes appears to be different 
and a flattening of the XLF of the radio-loud subsam- 
ple is visible for Lx{z = 0) < 3 x 10^** ergs^^. 

Ciliegi et al. (1995) have obtained VLA observa- 
tion of a ROSAT sample of 80 faint AGNs. They 
combined their sample with EMSS data and found 
that, in the case of a PLE, Lx{z) = Lx{0){l + z)'', 
fc = 2.43 ± 0.26 and k = 2.71 ± 0.10 for the radio-loud 
and radio-quiet population respectively. Although 
not significant, the difference hints a slower evolution 
of the radio loud population. 

In summary, at the moment data from different 
wavebands suggest that: 

1. There is a flattening of the bright part of the 
OLF at low redshift; 

2. There is an increase of the radio-loud fraction at 
Mb < —25 both in X-ray and optically selected 
samples; 

3. The apparent luminosity evolution of the radio- 
loud population is slower both in optically and 
in X-ray selected samples; 

4. The same QSO population is selected in the op- 
tical and soft-X domains, and X-ray selected 
QSOs follow a PLE evolution similar or slower 
than that of the optically selected QSOs; 



5. The X-ray selected samples do not have a suffi- 
cient coverage of the L, z plane to allow a mean- 
ingful cross-check of the flattening in the bright 
part of the low- 2: OLF. 

In this way the observed overdensity of low redshift 
bright QSOs and the related flattening of the LF could 
be originated by the contribution of a slower evolv- 
ing population of radio-loud QSOs. To put this sce- 
nario on a firmer basis it is necessary to obtain radio 
data for: 1) a new bright, wide field (~ 10000 deg^) 
QSO sample, in the interval 13 < B < 16, and 2) 
a wide field bright X-ray selected QSO sample. The 
former will be provided by the Sloan Survey (Gunn 
and Weinberg 1995), while the Rosat All Sky Survey 
seems to be ideally tailored for the latter. 
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TABLE 2. "Best fit" parameters for luminosity 
function models. 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 
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Table 1. QSO counts and surface densities. 



Magnitude range 
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Notes to Table 1. 

QSOs are defined as objects having Mb < —23 
and 0.3 < z < 2.2 with Hq = 50 Km/s/Mpc 
and qo = 0.5. 
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Table 2. "Best fit" parameters for luminosity function models. 



Model 


90 


mag 


$ 




fc2 
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P 


M* 


z' 


Prob. 
z < 2.2 


2DKS Prob. 
z < 0.6 
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-1.6 


-3.8 


-22.6 




0.004 


0.01 




0.1 


3.4 


X 10-^ 


3.44 




-1.57 


-3.75 


-22.9 




0.36 


0.08 


M ™ 


0.1 


4.0 


X lO-'^ 


3.45 


0.92 


-1.50 


-3.74 


-22.8 


0.40 


0.49 


0.11 


1 (7 errors: 








±0.07 


±0.27 


±0.07 


±0.13 


±0.2 









'^PLE model from Boyle (1992) which includes a Zcut — 1.9 parameter 

•^Our best fit PLE model 

=Our best LDLE model. 

'^PLE model without our HBQS data. 

°PLE model without the EQS (Goldschmidt et al. 1992) data. 

fpLE model without the LBQS (Hewett, Foltz and Chaffee 1995) data. 

SPLE model with a spectral slope dispersion — 0.3. 

''PLE model with a spectral slope dispersion a-y ~ 0.5. 

'PLE model from Boyle (1992) which includes a Zcut — 2.1 parameter 

'Our best fit PLE model 

•"Our best LDLE model. 
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